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Abstract 

We examine elastic Compton scattering off the deuteron for photon energies be- 
tween 50 MeV and 100 MeV in the framework of chiral effective field theories to next- 
to-leading order. We compare one theoretical scheme with only pions and nucleons as 
explicit degrees of freedom to another in which the A(1232) resonance is treated as 
an explicit degree of freedom. Whereas pion degrees of freedom suffice to describe the 
experimental data measured at about 70 MeV, the explicit A(1232) gives important 
contributions that help to reproduce the angular dependence at higher energies. The 
static isoscalar dipole polarizabilities a s E and (5 S M are fitted to the available data, giving 
results for the neutron polarizabilities cc^ = (14.2 ± 2.0 (stat) ± 1.9 (syst)) • 10~ 4 fm 3 , 
Pfo = (1.8±2.2 (stat) ±0.3 (syst)) -10" 4 fm 3 . These values are in good agreement with 
previous experimental analyses. Comparing them to the well-known proton values we 
conclude that there is currently no evidence for significant differences between the 
proton and neutron electromagnetic dipole polarizabilities. 
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1 Introduction 



The structure of protons and neutrons as analyzed with electromagnetic probes has been 
under experimental and theoretical investigation for a number of decades. In Compton scat- 
tering the external electromagnetic field of the photon attempts to deform the nucleon. The 
electromagnetic polarizabilities provide a measure of the global resistance of the nucleon's 
internal degrees of freedom against displacement in an external electric or magnetic field, 
which makes them an excellent tool to study the sub-nucleonic degrees of freedom. If one 
defines the polarizabilities via a multipole expansion in Compton scattering, these quantities 
are energy dependent. The physics connected with this energy dependence is discussed in 
Refs. 0121- In this work we are mainly concerned with the static values, which we therefore 
denote as the polarizabilities for simplicity. Experimentally the best-known nucleon polar- 
izabilities are the static electric and magnetic dipole polarizabilities of the proton, oP E and 
f3 p M . The results of the global fit to the wealth of Compton scattering data on the proton 
given in Ref. [B] are 

a p E = (12.1 ± 0.3 (stat) q= 0.4 (syst) ± 0.3 (model)) • 10" 4 fm 3 , 

(f M = (1.6 ± 0.4 (stat) ± 0.4 (syst) ± 0.4 (model)) • 10" 4 fm 3 . (1) 

The values from a Baldin Sum Rule constrained fit of the proton Compton data within the 
framework that we use in this work, reported in Ref. j2], agree within error bars with (JTJ: 

a p E = (11.04 ± 1.36) • 10" 4 fm 3 , (3 P M = (2.76 =f 1.36) • 10" 4 fm 3 . (2) 

The errors displayed in Eq. (J2J) are only statistical. For the fit, the central value of the 
Baldin Sum Rule a p E + f3 p M = (13.8 ± 0.4) ■ 10" 4 fm 3 [3] has been used. 

Due to the lack of stable single-neutron targets for Compton scattering it is much harder 
to access the neutron polarizabilities experimentally. An experiment on quasi-free Compton 
scattering from the proton and neutron bound in the deuteron [3] gives results for the 
neutron polarizabilities which suggest very small isovectorial components 1 when compared 
to Eqs. (P2): 

a n E = (12.5 ± 1.8 (stat) jjj; 1 (syst) ± 1.1 (model)) ■ 10" 4 fm 3 

f3^ = (2.7 =F 1.8 (stat) +°;° (syst) =f 1.1 (model)) ■ 10" 4 fm 3 (3) 

A similar observation for a E has been made in j^j , where the scattering of neutrons on lead 
was measured: 

a n E = (12.6 ±2.5) • 10~ 4 fm 3 (4) 

However, the precision of this result has been questioned by the authors of jH]. Their 
estimate of the correct range for the result from j^j is 7 < a' E < 19. On the other hand, 
another experiment jjj, using the same technique, gives a completely different result: 

a n E = (0.6 ± 5.0) ■ 10~ 4 fm 3 (5) 
1 The isovector polarizabilities are defined as a v E = \ (a p E — a E ), f3\ { = | (f3 p M — (3^)- 
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On the theory side Chiral Perturbation Theory predicts that the proton and neutron 
polarizabilities are equal at leading-one- loop order jH], since the pion loops that generate 
these contributions are isoscalar in nature. The absence of large isovector pieces in and 
13 m therefore is in accord with this picture. 

Another possible way to determine the neutron polarizabilities is elastic low-energy 
Compton scattering from light nuclei, e.g. from the deuteron. Several experiments have 
already been performed [§J EH E] and further proposals exist - e.g. Compton scatter- 
ing on the deuteron or He 3 at TUNL/HI7S jT^j and on deuteron targets at MAXlab [TB*] . 
From a theorist's point of view, extracting the neutron polarizabilities from elastic 7<i scat- 
tering requires an accurate description of the nucleon structure and the dynamics of the 
low-energy degrees of freedom within the deuteron, as one has to correct for the proton 
polarizabilities and meson-exchange effects. A first attempt to fit the isoscalar polariz- 
abilities a% = \ (a? E + a E ), Pm = \ [fi v M + 0m) to the elastic deuteron Compton scat- 
tering data from E] has been made in [Tl] . The extracted neutron polarizabilities 
a n E = (9.0 ± 3.0) • 10~ 4 fm 3 , pfo = (11.0 ± 3.0) • 10" 4 fm 3 indicate the possibility of a 
rather large isovector part. On the other hand comparison of the elastic deuteron Compton 
calculation of Ref. ^3] with the data from [5] is in good agreement with nearly vanishing 
isovector polarizabilities: a n E = (12.0 ± 4.0) • 10" 4 fm 3 , pfc = (2.0 ± 4.0) • 10" 4 fm 3 , albeit 
within rather large error bars. 

It is obvious from these partly contradictory results that there is still a lot of work 
to be done in order to have reliable values for a E and (3^. In general, Chiral Effective 
Field Theory provides a consistent, controlled framework for elastic 7c? scattering within 
which nucleon effects can be disentangled from meson-exchange currents, deuteron binding, 
etc. It therefore gives valuable contributions to the ongoing discussion on the neutron 
polarizabilities. This work aims for an improved description of elastic deuteron Compton 
data at uj ~ 50-100 MeV, compared to the calculations presented in [T4l I15j . which cannot 
match the data in this regime. 

Our work is based on the calculation of Refs. fHl EI, where Compton scattering off 
the deuteron was examined for photon energies u ranging from 50 MeV to 100 MeV. The 
central values for the isoscalar polarizabilities, derived in the recent analysis [Fj\ of the data 
from El E] are 

a s E = (8.9±1.5)l 4 j • 10~ 4 fm 3 , 

/^ = (2.2±1.5)+^10- 4 fm 3 . (6) 

Comparing with Eq. J!}, these results indicate a small isovector magnetic polarizability, 
but signal the possibility of a rather large ot v E . However, the range for a s E and (3 S M quoted 
in PH is rather large: a s E = (8.0 - 13.6) • 10~ 4 fm 3 , (3 S M = (1.3 - 3.4) • 10~ 4 fm 3 . The 
authors of Ref. ^H] followed Weinberg's proposal to calculate the irreducible kernel for the 
ryNN — > 'yNN process in Heavy Baryon Chiral Perturbation Theory (HB^PT) and then fold 
this with external deuteron wave functions such as Nijm93 , CD-Bonn ^1 or AV18 [2*Uj . 
Proceeding in this fashion means working within an Effective Field Theory in which only 
nucleons and pions are active degrees of freedom. This "hybrid" approach has proven quite 
successful in describing ird e ~d |22j . and even jd ^21 Ej scattering. In this work we 
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extend the calculation of Ref. ^H] to an Effective Field Theory which includes the A (1232) 
resonance of the nucleon as an additional explicit degree of freedom. The advantage of our 
approach with respect to the NNLO calculation of Ref. is that we identify the physics 
hidden in some of the short-distance parameters there. As we shall see, this is particularly 
important for quantities such as the magnetic polarizability (3m, where the A(1232) plays 
an important dynamical role. The huge influence of the A (1232) in single- nucleon Compton 
scattering - especially in the backward direction - is a well-known fact (see e.g. [2]). We note 
that the importance of the A-contributions is due to the strong paramagnetic Ml coupling 
of the photon to the N — > A transition, visible already far below resonance (cf. j2]). It is 
therefore interesting to also investigate the role of these degrees of freedom in elastic 'yd 
scattering, and this is the main focus of our work. 

The idea of an extension of HB^PT that includes explicit A degrees of freedom has its 
origin in the early 1990's [22]. When including the A(1232) explicitly in xEFT one needs 
to specify how the AN mass splitting Ao is treated in the power counting. Here we use 
the so-called Small Scale Expansion (SSE) [21]. We note that there also exist alternative 
approaches for Chiral Effective Field Theories with explicit tt, N and A degrees of freedom, 
e.g. the ^-expansion [53], which was recently shown to describe jp cross section data well 
in an energy range from uj = MeV to u — 300 MeV. 

In Sect. E] we discuss our predictions for the deuteron Compton cross sections for three 
different energies between 50 MeV and 100 MeV, comparing to experimental data and to 
the 0(q 3 ) HB^PT calculation ^H]. Before that, we give a brief survey of the theoretical 
formalism in Sect. [2] There we show that combining Weinberg's counting ideas with the 
SSE power counting scheme leads to no additional diagrams in the two-body part of the 
kernel with respect to [IB]. In Sect. [U we present our results for the isoscalar polarizabilities, 
derived from a fit to elastic deuteron Compton scattering data, which turn out to be in 
good agreement with the theoretical expectation that the isovector components are small. 
We conclude in Sect. El and give a brief outlook on future projects, one of which aims to 
cure the shortcomings of our calculation in the extreme low-energy regime uj <C 50 MeV (cf. 
Sect. 12]). 

2 Compton Scattering off the Deuteron in Effective 
Field Theory 

We are calculating Compton scattering off the deuteron in the framework of the Small Scale 
Expansion [21], an Effective Field Theory with nucleons, pions and the A(1232) resonance 
as explicit degrees of freedom. In this extension of xPT, the expansion parameter is called 
e, denoting either a small momentum, the pion mass or the mass difference A between 
the real part of the A mass and the nucleon mass. The relevant pieces of the Ntt and An 
Lagrangeans have been discussed in the literature many times, and we refer the interested 
reader to [2E| for the Nn Lagrangean and for a general review of HB^PT, and to [21] and [2] 
for the relevant pieces of the An Lagrangean. 

The power-counting scheme that we use for Compton scattering off light nuclei is mo- 
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tivated by Weinberg's idea to count powers only in the interaction kernel. We base our 
calculation on the hybrid approach, which is a well-established tool by now. While the 
kernel is power counted according to the rules of the Effective Field Theory, the deuteron 
wave functions we use are obtained from state-of-the-art NN potentials: Nijm93 [IB], the 
CD-Bonn potential [IE], the AV18 potential [20] and the NNLO chiral potential [2B], where 
this last potential also follows Weinberg's suggestion, and is derived by applying HB%PT 
power counting to the NN potential V. 

It is convenient to write the Green's function for Compton scattering from the two- 
nucleon system as 

G?77 = G .fCy-y G -\- G Kry G G , (7) 

with i^ 77 the two-nucleon-irreducible part of the interaction kernel, which contains both one- 
body and two-body mechanisms, and G K 1 the so-called nuclear-resonance contribution 2 
to the kernel. G is the two-particle Green's function, constructed from the two-nucleon- 
irreducible interaction V and the free two-nucleon Green's function. We apply the same 
power counting rules to both K~ n and K^GK^, calculating all contributions up to a specific 
order in e, namely e 3 . 

To do this we first note that a diagram contributing at a certain order in q in HB^PT 
contributes at the same order e in SSE. However, there are some kinematics in which HB^PT 
counting should not be employed for the ^NN — > '-/NN kernel. In HB^PT the leading 
order propagator of a nucleon with the energy uj of the external probe flowing through it 
is — . Corrections from the kinetic energy of the nucleon are treated perturbatively. In the 
deuteron, such a perturbative treatment is not applicable for low photon energies, due to the 
relative momentum p of the two nucleons. Therefore, in the low-energy regime one has to 
use the full non-relativistic nucleon propagator u _ p y 2 M • Nonetheless, the approximation ^ 
is useful for uj 3> B, with B m 2.225 MeV the binding energy of the deuteron, as p 2 /M ~ B 
with a typical nucleon momentum p inside the deuteron. These considerations demonstrate 
that for uj ^> B the nucleon propagator may be counted as 0{e~ 1 ) like in standard HB^PT, 
whereas in the 'nuclear' regime uj ~ 0(m 2 /M) it has to be counted as 0(e~ 2 ), as p ~ e. 
Therefore, one has to strictly differentiate between two energy regimes: the nuclear regime 
uj ~ 0(m 2 /M) and the regime uj ~ 0(m n ). Here we work in the latter one, as we are 
mainly concerned with photon energies uj > 50 MeV, which is the energy region where one 
starts to be sensitive to the nucleon polarizabilities. 

We note further that only in the regime uj ~ 0(m n ) can one treat the contributions from 
K^GK^ using a perturbative chiral expansion. Since we do treat this piece using HB%PT 
counting it is therefore no surprise that the Compton low-energy theorems are violated. For 
example, the Thomson limit for Compton scattering from a nucleus of charge Z e and mass 
AM, 

Z 2 e 2 _ 

^Thomson = A((J = 0) = -— 6 ■ 6* , (8) 

a direct consequence of gauge invariance [2H] , cannot be recovered without the full resonance 

2 The nomenclature is due to nuclear resonances which are excited by the initial interaction with the 
photon and which one might expect to dominate the Compton process at low energies, see e.g. |14) . 
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term. Therefore, we strictly constrain ourselves to photon energies u ~ 0{m 7T ), where a 
perturbative expansion of the kernel in the standard HB^PT counting scheme, i.e. counting 
the nucleon propagator as (9(e _1 ), is possible. The lower limit of this power counting turns 
out to be uj ~ 50 MeV, so we have to caution the reader that the calculation is not supposed 
to work in the region uj <C 50 MeV (an extension to lower energies is in progress |SD])- A 
more detailed discussion of the power counting applied to the meson exchange part of our 
calculation can be found in Refs. [To^ li7|. For an Effective Field Theory approach to deuteron 
Compton scattering where pions are integrated out, see Ref. jSIj. These calculations describe 
the very- low-energy region well and also reach the exact Thomson limit. 

The T-matrix for Compton scattering off the deuteron is derived as the matrix element of 
the interaction kernel, evaluated between an initial and final state deuteron wave function, 
as explained in great detail in [To] . 



Stated differently, one obtains T by extracting the piece of G 77 corresponding to the deuteron 
pole at E = —B in both the initial and final state. 

As we are calculating 'yd scattering in the Small Scale Expansion, we also have to fix our 
counting rules for diagrams including A(1232) propagators. For the one-body contributions 
this is straightforward, as we apply the SSE counting scheme, cf. Refs. [21 E2]- As far as 
the two-body physics is concerned, we combine the SSE counting rules, e.g. counting the 
A-propagator as e -1 , with Weinberg's prescription of counting only within the interaction 
kernel. To 0(e 3 ), the order up to which we are working, this leads to identical meson 
exchange diagrams as in the 0(q 3 ) HBxPT calculation. All additional diagrams are at 
least one order higher, an example is given in Fig. E^b) (an example of an 0(e 4 ) one-body 
diagram is sketched in Fig. EJa)). A modified counting scheme in the two-body sector has 
been suggested in j2I] , as certain pion-exchange diagrams may be enhanced when the photon 
energy comes close enough to the pion mass that the pions in the two-body diagrams are 
almost on mass shell. We do not consider such a modification necessary for our calculation, 
as we restrict ourselves to photon energies uj < 100 MeV. 



Figure 1: Two examples of 0(e ) contributions to Compton scattering on the deuteron 
including explicit A(1232) degrees of freedom for one- and two-body contributions. 

Therefore, the diagrams contributing to deuteron Compton scattering up to 0(e 3 ) are: 

• One-body contributions without explicit A(1232) degrees of freedom. These are the 
single-nucleon seagull with the two-photon vertex from (Fig. Efa)), which is the 



T= (VflK^ + K^GK^Vi) . 



(9) 




(a) 



(b) 
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only 0(e 2 ) contribution, the nucleon-pole terms (Fig. 12(b)) , the pion pole (Fig. 12(c)) 
and the contributions from the leading chiral dynamics of the pion cloud around the 
nucleon (Figs. |2(d)-(g)). Up to third order, the only difference in these diagrams, 
compared to Compton scattering off the single nucleon, is that the pole diagrams 
(Fig. E(b)), which are conveniently calculated in the 77V center of mass frame, have 
to be boosted to the 'yNN center of mass system, as the calculation is performed in 
the 'yd cm frame. The resulting formulae from the boost are given in |16j . 
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(d) 



(f) 



(9) 



Figure 2: One-body interactions without a A(1232) propagator contributing to Compton 
scattering on the deuteron up to 0(e 3 ) in SSE. Permutations and crossed graphs are not 
shown. 



One-body diagrams with explicit A degrees of freedom, as shown in Fig. El The A-pole 
diagrams (Fig. 0(a)) and the contributions from the pion cloud around the A(1232) 
(Figs.O(b)-(e)). 

Two isoscalar short-distance one-body operators (Fig. 0(f)), which give energy-inde- 
pendent contributions to the dipole polarizabilities a s E and (3 S M . They are formally of 
(9(e 4 ) but turn out to give an anomalously large contribution to the single-nucleon 
Compton amplitude. Therefore, they have to be promoted to next-to-leading order 
as discussed in detail in j2], which we also refer to for the Lagrangean. We note that 
except for the two contact operators (Fig. 0(f)), the 5-expansion (cf. Sect. [TJ up to 
NNLO is equivalent to 0(e 3 ) SSE in the energy range ui ~ considered. 



(a) 



(b) 




V, 




(d) 



(e) 



(f) 



Figure 3: One-body interactions which contribute to deuteron Compton scattering at 0(e 3 ) 
in SSE in addition compared to third-order HBxPT. Permutations and crossed graphs are 
not shown. 
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Two-body contributions with one pion exchanged between the two nucleons (Fig. EJ). 
In total there are nine two-body diagrams at 0(e 3 ). As discussed before, the meson 
exchange diagrams are identical in third-order HB^PT and SSE. 



(a) 



(e) 



Figure 4: Two-body interactions contributing to the kernel for Compton scattering on the 
deuteron at 0(e 3 ) in SSE. Diagrams which differ only by nucleon interchange are not shown. 



All these diagrams (Figs. EH3| make up our interaction kernel. The SSE single-nucleon 
amplitudes can be found in [2j, while the two-body contributions are given explicitly in [Tfij . 
Note that we have simplified the expressions given in [2j with respect to the exact position 
of the pion threshold as we are only analysing Compton scattering for photon energies 
< 100 MeV. An estimate of the (small) size of this simplification is given in Sect. 13.21 

In the next section we compare our 0(e 3 ) SSE results for the deuteron Compton cross 
sections to the 0(q 3 ) HBxPT calculation performed in ^H] and to the available experimental 
data. Special interest is put on the energy and wave-function dependence of the cross section. 



3 Predictions for Deuteron Compton Cross Sections 

In Fig. El we compare the C(e 3 ) SSE predictions to the 0(q 3 ) HB%PT calculation of Ref. [Ej, 
using the wave function derived from the NNLO chiral potential with spectral function regu- 
larization [2H]- We also show the 0(q 2 ) result, which consists only of the single-nucleon seag- 
ull (Fig. 12(a))- The experiments shown have been performed at a lab-energy of 49 MeV [5], 
~67 MeV PII|, 69 MeV and ~94.2 MeV [HI- (The last experiment used photons in an 
energy range from 84.2 — 104.5 MeV; the deviation from the central value has been corrected 
for PI].) 

The values we use for physical constants can be found in Table ^ For the coupling 
constants connected with the two short-distance 7iV-operators (cf. Sect. EJ) and the jNA 
coupling bi we use the results of the Baldin Sum Rule constrained fit to the spin-averaged 
proton Compton scattering data from [2|. Fitting the short-distance operators is equivalent 
to fitting the static polarizabilities a E , f3 p M . We use the central values of the fit, which are 
a p E = 11.04 • 10" 4 fm 3 , (3 P M = 2.76 ■ 10" 4 fm 3 0, cf. Sect. □ The C(e 3 ) SSE calculation 
then predicts a E = a E , (3^ = [3 P M as the isovector contributions only come in at C(e 4 ). 
Therefore, there are no free parameters in our deuteron Compton calculation. 

From the 49 and 69 MeV curves shown in Fig. El it is obvious that explicit A degrees of 
freedom may well be neglected for these low energies. The two calculations - HB^PT and 
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Figure 5: Comparison of the 0(q 3 ) HB^PT (dashed) and the 0(e 3 ) SSE (solid) prediction 
at a;i ab = 49 MeV, u }ah = 69 MeV and w lab = 94.2 MeV using the chiral NNLO wave 
function [2E|. The data are from Illinois [H] (circle), Lund [TU] (star) and SAL [TTJ (diamond). 
The dotted line is the 0(q 2 ) result. 



SSE - yield results which differ only within the uncertainties one expects from higher order 
contributions. This is an important check, as it demonstrates the correct decoupling of the 
resonance, which provides the same low-energy limit in both theories. 

Another interesting observation can be made: The counting scheme described in Sect. |2] 
seems to break down for energies somewhere near 50 MeV, as both theoretical descriptions 
miss the 49 MeV data points, whereas the 69 MeV data (we neglect the minor corrections 
due to the data of [TU] being measured around 67 MeV) are well described within both 
theories. The 49 MeV data are best described by the 0(q 2 ) calculation but we believe 
this is a coincidence, as the low-energy theorems are violated at this order too. However, 
for higher energies - i.e. for describing the 94.2 MeV data correctly - the inclusion of the 
explicit A field seems to be advantageous in a third-order calculation. Here, 0(q 3 ) xPT 
misses the data in the backward direction. It also fails to reproduce the shape of the data 
points, which shows a slight tendency towards higher cross sections in the backward than in 
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Parameter 


Value 


Comment 


m n 


loy.o Mev 


charged pion mass 


M 


no o n i\/r i r 

938.9 Mev 


isoscalar nucleoli mass 


Jtt 


QO /I T\/ToV 

yz.4 ivie v 


pion decay constant 


9a 


1 1^7 


axial coupling constant 


Ot 




Ill-ill t~i no cf T'li r>T 1 1 VQ r~> r\ n o i o n i 

Vc^H/Ly llllt! SLIUClUie COllbLclllL 


K v 


o. (Uo 


isovector anomalous magnetic moment 


K s 


n us 

— U. 1 lo 


isoscaiai anomalous magnetic moment 


Md 


lo ( o.oo ivie v 


deuteron mass 


r> 
D 


Z.ZZ4D Mev 


deuteron binding energy 


A 

^0 


971 1 MoV 

zi i.i ivie v 


iV^A mass spiiiimg 


gnNA 


1.125 


7riVA coupling constant 


h 


4.67 


7iVA coupling constant 


gin 


18.82 


short-distance coupling constant 


0118 


-6.05 


short-distance coupling constant 



Table 1: xEFT parameters determined independently of deuteron Compton scattering. 
Magnetic moments are given in nuclear magnetons. 

the forward direction. This shape is very well reproduced in SSE, demonstrating once again 
the importance of the A resonance in Compton backscattering, due to the strong Ml —> Ml 
transition. We note that this feature can be clearly seen in the dynamical magnetic dipole 
polar izability (3mi{w), even for photon energies below the pion-production threshold (cf. [2] 
and Sect. 14. Calculations like the ones presented in Refs. [T4"l UK] , which include the 
dynamics of the polarizabilities only via the leading f3J and subleading terms [H] of a 
Taylor expansion, may therefore fail to describe the data around 95 MeV. 

3.1 Energy Dependence of the Cross Sections 

In order to decrease the statistical uncertainties, the experiment had to accept scattering 
events in an energy range of around 20 MeV. Therefore we think it worthwhile to examine 
the sensitivity of our results to the photon energy. In fact, our calculations suggest that 
the forward-angle cross section, in particular, has a sizeable energy dependence, which is, 
however, nearly linear. In Fig. |U] we show our results for three different photon energies 
around 69 MeV and 94.2 MeV, respectively, separated only by 10 MeV. This emphasizes 
the importance of having a well-defined photon energy at which to examine the effects of 
cce and /3m- 

3.2 Correction due to the Pion-Production Threshold 

In low-order HBxPT/SSE calculations the 7<i — > nNN threshold is not at the correct 
position as dictated by relativistic kinematics. For a similar problem, regarding the correct 
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0U [deg] #iab [deg] 

Figure 6: C(e 3 ) SSE results for 64 MeV, 69 MeV, 74 MeV and, respectively, 90 MeV, 
94.2 MeV, 100 MeV (from the upper to the lower curve in each panel), using the xPT wave 
function |2~H] . 




3 5 
2.5 



25 50 75 100 125 150 175 
^ab [deg] 

Figure 7: Estimate of the effect of a threshold correction (dotdashed) on the C(e 3 ) SSE 
results (solid), using the chiral NNLO wave function Im- 



position of the pion-pro duct ion threshold in the single-nucleon sector, see e.g. 0. So far we 
refrain from an analogous correction for scattering. However, in Fig. [7| we investigate 
what deviations one would expect from our present results, as indicated by an estimate, 
where we use the single-nucleon SSE amplitudes j2] with the exact expression for ^fs — M. 
Obviously, even at the highest photon energies considered here, 94.2 MeV, the corrections 
are negligible, given the sizeable error bars of the experimental data and the theoretical 
uncertainties of a leading-one-loop order calculation. 

3.3 Wave-Function Dependence of the Cross Sections 

Another interesting issue is the wave- function dependence of our results. Fig. |H] investigates 
the sensitivity to the wave function chosen, showing sizeable deviations between the NNLO 
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0(e 3 ) SSE results for 69 and 94.2 MeV, using four different wave functions: NNLO 
(solid), Nijm93 [HI (dotted), CD-Bonn [H] (dashed), AV18 [2U] (dotdashed). 



xPT wave function [2E] on one hand and the wave functions derived from the Nijm93 
potential [TH] and the AV18 potential [20] on the other. The last two yield results which 
are nearly indistinguishable but are considerably higher than those found with the wave 
function of (2Hj- With the CD-Bonn wave function JH] we obtain results in between NNLO 
%PT and Nijm93/AV18. This pattern is identical for both energies under investigation, 
69 MeV and 94.2 MeV. Given that our calculation is based on a low-energy Effective Field 
Theory of QCD, the dependence on the wave function is somewhat worrisome and will 
be discussed further in jSU]. According to Weinberg counting it is an 0(q 5 ) effect, so a 
deviation of the order of 10% is more than one would expect. We interpret this feature 
as an unwanted sensitivity to short- distance physics, because the long-range part of all 
wave functions, described by one- and two-pion exchange, is identical. However, one must 
caution that the NNLO xPT potential reproduces the Nijmegen partial-wave analysis with 
less precision than the CD-Bonn, AV18 or Nijm93 potentials. 

In this section we presented our predictions for differential cross sections. These are 
parameter-free as we fixed the nucleon polarizabilities via proton Compton data. The good 
agreement of the SSE results with experiment at 69 MeV and 94.2 MeV leaves little room for 
large isovector polarizabilities, since these predictions used the same values for the proton 
and neutron polarizabilities. It further encourages us to determine the isoscalar dipole 
polarizabilities a s E and (3 S M directly from the deuteron Compton cross sections. The results 
are displayed in the next section, together with the results one obtains from analogous fits 
using the 0(q 3 ) HB^PT amplitudes. 



4 Determining a s E and f3 s M from Scattering 

An accurate and systematically-improvable description of Compton scattering on deuterium 
offers the possibility to extract the isoscalar polarizabilities directly from deuteron Compton 
scattering experiments in a systematic way. The resulting numbers can then be combined 
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with the known numbers for the proton to draw conclusions about isovector pieces a v E and 
f3 v M , or, equivalently, the elusive neutron polarizabilities. As our SSE calculation provides 
a reasonable description of the 69 MeV and the 94.2 MeV data (see Sect. |3J), we present in 
the following our results from a fit of the isoscalar polarizabilities to these two data sets. 
This corresponds to fitting the coupling strengths of the two short- distance isoscalar 7iV- 
operators (Fig. Eff)), which we now fit to 'yd data rather than to 7p data. In this way we 
can check our assumption that the short- distance operators are isoscalar at leading order. If 
the value extracted from 'yd data is approximately that from 7p data, that argues in favour 
of short- distance mechanisms which are predominantly isoscalar. 

Our SSE results are compared to the fit results that we get for a s E and j3 s M when we use 
modified 0(q 3 ) HBxPT amplitudes. This modification consists of including in our calcu- 
lation isoscalar short-distance jN operators which change both the electric and magnetic 
polarizability from their 0(q 3 ) values. In other words, we write 



5 a 9 a 

96 /2 mir tt 



* = i927fb +, <'- (10) 

The energy dependence of the polarizabilities is still given solely by the leading-order pion 
cloud. Eq. (JTUJl promotes the short-distance contribution to a and f3 from 0(q 4 ) to 0(q 3 ). 
There are indications that this change in the power counting is necessary if high-energy 
modes in the pion- loop graphs that generate a and j3 are to be properly accounted for 
In order to avoid confusion we denote the fits done with this procedure as HBxPT 0(q 3 ). 

Fits similar to our 0(e 3 ) and 0(q 3 ) ones have already been performed in [TJJ, calculating 
in HB^PT up to 0(q A ). The authors of [T7] used all available data sets but had to exclude 
the two 94.2 MeV data points measured in the backward direction. As jTHl EI and the 
0(e 3 ) SSE calculation obviously have problems to describe the normalization of the data 
at energies below 60 MeV we decided to only include the data around 69 MeV EH] and 
94.2 MeV in the fit. We do not make any cuts on the angles and, in contradistinction 
to ^7], we do not allow the normalizations in the various experiments to float in the fit 
within their quoted systematic errors. 

We performed the fits using the NNLO chiral wave function. We fitted the 16 data 
points using 2 free parameters (a s E and f3 s M ), leaving us with 14 degrees of freedom. The 
resulting values for a% and (3 S M (see Table |2J) are 

a s E = (12.1 ± 1.3) • 10" 4 fm 3 , (3 S M = (1.8 ± 1.6) • 10" 4 fm 3 (11) 

with a x 2 /d.o.f. of 1.78. The corresponding plots are displayed in Fig. El together with 
the results of our 0(q 3 ) fits. Using the experimental values from Eq. JEJ [3] as input 
or, equivalently, the values given in Eq. J2J), which are obtained from proton Compton 
scattering data within the same framework as we are using here, one can derive the neutron 
polarizabilities from the isoscalar ones: 

a n E = (12.1 ± 1.3) ■ 10" 4 fm 3 , = (2.0 ± 1.6) • 10" 4 fm 3 (12) 
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From these results we deduce that the isovector polarizabilities are rather small (see Table[2J), 
in good agreement with expectations, which predict the isovector part to be of higher 
than third order. Therefore we find no contradiction between the results from quasi-free [3] 
and elastic deuteron Compton scattering. 

Our results for a s E and f3 s M in SSE (cf. Eq. (jllj) ) are well consistent (within error bars) 
with the isoscalar Baldin Sum Rule 



a s E + (3 s 



M 



= (14.5 ±0.6) • 10" 4 fm 3 , (13) 

world av. 



which has been a serious problem in former extractions [T4l I17j . The numerical value for 
the sum rule is derived from 

a p E + f3 p M = (13.8 ± 0.4) • 1(T 4 fm 3 [3], 

a% + 0h= (15.2 ±0.5)- 10" 4 fm 3 @|. (14) 

Due to the consistency of our fit results with the sum rule value from Eq. (fT3*|) one can in a 
second step use this number - we use the central value - as an additional fit constraint and 
thus reduce the number of free parameters to one. The resulting one-parameter fits in SSE 
of Table H 

a% = (12.3 ± 0.7) • 10" 4 fm 3 , (3 S M = (2.2 t 0.7) ■ 10~ 4 fm 3 , (15) 

are in good agreement with the isoscalar average of the numbers from Eqs. (0) and (jSJ) - or, 
alternatively, Eqs. (JTJ) and (HJ). 

Comparing our fit results to the isoscalar 0(q 4 ) HB%PT estimate [33] , &% = (H-95 ± 
2.5) • 10~ 4 fm 3 , p s M = (5.65 ± 5.1) • 10~ 4 fm 3 , we see only minor deviations from their value 
for a s E . However, our values for f3 M are significantly smaller, but no meaningful conclusion 
can be drawn due to the large error bars in the 0(q A ) estimate. The reason for the huge 
error bars in the 0(q A ) HBxPT numbers of Ref. [33] is their sensitivity to short-distance 
contributions which were estimated using the resonance-saturation hypothesis. 



4.1 Wave-Function Dependence of the Fits 

To have an estimate on the systematic error due to the wave-function dependence, we show 
our results when we use the two extreme wave functions (cf. Fig. |HJ) for the fit: the NNLO 
chiral wave function [2H] and the wave function from the Nijm93 potential [TSj. Furthermore, 
we are fitting in two different ways: First the number of degrees of freedom is the number 
of data points (16) minus the number of free parameters (2). In a second step we use the 
isoscalar Baldin Sum Rule, Eq. (JT3J), to reduce the number of degrees of freedom to 1, as 
described before. 

Fitting the 'yd cross sections with the 0(e 3 ) and 0(q 3 ) kernel, respectively, using the 
Nijm93 wave function yields larger results for a# and smaller ones for (3m, but still the values 
of Table El are in reasonable agreement with the values given in Eq. (J3j) [3] . Comparing the 
differing results that we get for a s E with the Nijm93 and the NNLO wave function, we 
estimate our systematic error to be of the order of 15 %. 
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Table 2: Values for the isoscalar and neutron polarizabilities from a fit to the full 69 MeV 
and 94.2 MeV data sets 011011111, using the C(e 3 ) SSE amplitudes and the 0(q 3 ) HByPT 
amplitudes, respectively. The neutron results are derived using the proton values from [3] 
(Eq. (UJ) as input. All error bars displayed are only statistical. 



One of the reasons for the differing results between our approach and the calculations 
presented in J2J is the energy dependence of the polarizabilities. In our calculation, it 
is completely given by the Compton multipoles, whereas the authors of f3J only used 
the leading [TU [T3] and subleading [T3] terms of a Taylor expansion of 0^(0;) and Pm(uj) 
in the photon energy uj. Comparing the range of our results for a s E and (3 S M to the ranges 
quoted in [T7] (cf. Sect. [TJ), we observe a slight tendency towards larger values of a s E and 
smaller ones of (3 S M in our analysis. The reason for this deviation will be discussed in detail 
in Sect. 0~1 

4.2 Comparison of C(e 3 ) SSE and 0{f) HBxPT Fits 

When we compare the 0(e 3 ) SSE fit results for a s E and j3 s M with the corresponding 0(q 3 ) 
HByPT results (Table 0), we see that in the HByPT fit (Eq. (JTj|) the electric dipole 
polarizability is smaller, whereas (3 S M turns out to be larger. The reason for the systematic 
shift of the magnetic polarizability is that due to the missing A(1232) resonance in HByPT 
the static value of (3m is inflated in order to compensate for the paramagnetic rise of the 
resonance. This will be discussed further in the next subsection. 

Fig. El demonstrates that this compensation works very well in the 70? cross sections, as 
the curves, which correspond to the 0(e 3 ) SSE and to the 0(q 3 ) HByPT fits, are nearly 
indistinguishable. We consider the quality of our 0(q 3 ) fit to be comparable to that of the 
0(q 4 ) fits of Ref. [J7]. 

The SSE and HByPT fit results shown in Fig.|2]only differ in the associated pairs a^, (3m- 



15 



25 




<x) lab =69MeV 


B 20 






J 15 
G 




i N — ^PTt 


b 10 




T3 




W 5 





25 50 75 100 125 150 175 
flab [deg] 



25 



B 20 

J 15 
G 

b 10 

^ 5 





Wlab =69MeV 






1 






25 50 75 100 125 150 175 
ftab [deg] 




25 50 75 100 125 150 175 
0i ab [deg] 



25 50 75 100 125 150 175 
flab [deg] 



Figure 9: £>(e 3 ) SSE (solid) and 0{f) HB^PT (dashed) results with oe%, (3 S M from Table H 
using the chiral NNLO wave function [2E|. The upper panels correspond to a fit of both 
polarizabilities, in the lower panels the Baldin Sum Rule (cf. Eq. (JIEJ)) is used as additional 
fit constraint. The grey bands are derived from our (statistical) errors. 



Therefore, from the available 70? data alone one cannot draw any firm conclusion regarding 
the importance of explicit A (1232) degrees of freedom. However, as we will show in the 
next section, from 7^ scattering experiments it is clear that third-order HB^PT does not 
describe the dynamics in the 7c? process correctly. Given that the SSE calculation describes 
both the 7p and the 7c? experiments we believe we have established that a Chiral Effective 
Field Theory which includes the explicit A field is an efficient framework with which to 
describe low-energy Compton scattering. 

4.3 Towards an Unbiased Fitting Procedure 

When we compare the several data sets that we use for the fits, we find eleven data points 
at u;i a b ~ 69 MeV, centered around only two different angles, and five points at ~ 
94.2 MeV, distributed over the whole angular spectrum. Especially around # lab ps 130° 
there is a wealth of data around 69 MeV (six points from ^U] and one from [H]), which 
gives an anomalously strong constraint to our fit routines. As long as there are no further 
data available at higher energies, fitting to all of the 69 MeV data thus overestimates this 
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angle [deg] 


da/dQ [nb] 


45.6 
130.5 


17.3 ±2.8 
14.2 ±2.1 



Table 3: Effective data points representing the 67 MeV data from |10j . 



energy region with respect to the 94.2 MeV data from [TT]. Therefore, in the following we 
compensate for this imbalance of data by replacing the Lund data ^U] by two "effective" 
pseudo-data points (cf. Table El and Fig. ITTj) . which represent the data in the forward 
and backward direction, respectively. These are obtained by weighting the angles and the 
differential cross section values of the represented data points by the inverse of their errors, 
and we assign the average over the errors of the represented data as error bars. Therefore, 
the remaining data are the two data points from jU] at 69 MeV, the two "effective" data at 
~67 MeV, shown in Table El representing ^U], and the five data points from around 
94.2 MeV. With these nine data points we perform the same fits as we did before for the 
complete data sets. The resulting values for a# and (3 m are presented in Table El the plots 
(including the two effective data points) in Fig. El exhibiting better agreement with the 
94.2 MeV data than the fits of Sect. 14.21 (Fig. as expected. Comparing the results for a s E 
and (3 S M of Table |U (or, equivalently, for a E and (3 M ) to the results from our fits to all data 
points, given in Table El we note that both for 0(e 3 ) SSE and 0(q 3 ) HB^PT o>e is slightly 
smaller {(3m slightly larger). 

We also see once again that the theory without explicit A degrees of freedom leads to a 
systematically larger value for (3 M , supporting our hypothesis that the enhancement is due to 
the insufficient dynamics in the HBxPT (cf. Tabled- This is demonstrated in Fig. EH which 
shows the dynamical isoscalar dipole polarizabilities a E1 (u), (3' M1 (u), calculated from third- 
order SSE and HB^PT, respectively. The static values are taken from the unconstrained 
fit, using the NNLO %PT wave function, i.e. 

a% = 11.5-10- 4 fm 3 , (3 S M = 2.4 • 10" 4 fm 3 (16) 

for the C(e 3 ) SSE-curve, and 

a s E = 10.4 • 10~ 4 fm 3 , (3 S M = 3.5 ■ 10~ 4 fm 3 (17) 

for the 0(q 3 ) %PT -curve. As shown in Fig. [HI] the energy dependence of the two field- 
theoretical calculations for cxei{oj) is in good agreement with each other and with the recent 
analysis from Dispersion Theory [2|. Matters are different for (3mi{uj): Whereas the SSE- 
curve reproduces the paramagnetic rise due to the explicitly included A resonance, the 
HB^PT result amounts to a nearly energy-independent average for this quantity. However, 
it is well-known that this paramagnetic rise in the Compton multipoles is necessary for 
the correct description of the 7p data around pion threshold, as can also be seen in the 
Dispersion- Theory analysis. In HBxPT the static value is artificially enhanced by the fit 
constraint from the 94.2 MeV data in order to compensate for the missing dynamics. Both 
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in Fig. Eland Fig. ITDthe enhanced (5m is able to cure the 70? cross sections and make the 
resulting curves very similar to the plots from the SSE fits (cf. Fig. fTT|) . 

Therefore, for understanding the available 7c? data via fits of a s E and (5 M , it is essential to 
combine the pairs ag, (5m resulting from the 7<i analysis with an energy-dependent multipole 
analysis of scattering. From the information available on Compton multipoles from 7^ 
scattering experiments, it is clear that third-order HBxPT is too simplistic a picture for the 
dynamics of the 7c? process at energies of (9(100 MeV). It is therefore crucial that deuteron 
and proton Compton experiments are available at comparable energies and that they are 
analyzed within the same framework. 

Putting equal statistical weight on the 69 and the 94.2 MeV data can be seen as a demon- 
stration of the importance of obtaining comparable statistics at all energies. We therefore 
urge for more experimental information at photon energies around 100 MeV. With such 
information, deuteron Compton cross sections below the pion mass provide an excellent 
window to investigate which internal nucleonic degrees of freedom contribute in both pro- 
cesses, 7p — > 7p and 7c? — > 7c?. 




20 40 60 80 100 120 140 20 40 60 80 100 120 140 

w [MeV] w [MeV] 

Figure 10: Comparison of the isoscalar Dispersion Theory result (dotted) of Ref. [2| for the 
dynamical dipole polarizabilities a s E1 (uj) and (5 s mi {uj) to the C(e 3 ) SSE (solid) and the 0(<f) 
HB^PT (dashed) results with the static values from the two-parameter fit using the chiral 
wave function, given in Table EJ uj^ denotes the pion-production threshold. 



5 Conclusion and Outlook 

In this work Compton scattering from the deuteron was calculated up to next-to-leading 
order in the Small Scale Expansion, an Effective Field Theory with nucleons, pions and the 
A(1232) resonance as explicit degrees of freedom. We investigated three different energies 
and compared to the available experimental data, finding good agreement for 69 MeV and 
94.2 MeV and a failure of our calculation for 49 MeV. The reason for this last result is that for 
low energies our power counting breaks down, so one has to modify the power counting in the 
very-low-energy region and use non-perturbative methods in order to reproduce the correct 
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Table 4: Values for the isoscalar and neutron polarizabilities from a fit to the 69 MeV 
and 94.2 MeV data sets [TT], using the C(e 3 ) SSE amplitudes and the O(f) HByPT 
amplitudes, respectively. The data from [TU] have been replaced by two pseudo-data points, 
specified in Table El The neutron results are derived using the proton values from [3] as 
input. All error bars displayed are only statistical. 



low-energy theorems. It will be one of our future projects [3U] to get the region uj < 50 MeV 
under control and finally to restore the correct Thomson limit. Here we concentrated on the 
energy range between 50 MeV and 100 MeV. We found that our calculation gives reasonable 
results for photon energies above 60 MeV. Motivated by the good agreement at these higher 
energies, we fitted the isoscalar polarizabilities a s E and (5 S M to the data around 69 MeV and 
94.2 MeV, yielding results in good agreement with the 0(q 4 ) HByPT estimate of Ref. 
and experiment. Averaging over the results of our two unconstrained SSE fits (one with 
the chiral NNLO wave function j2H], one with the Nijm93 wave function [T^j, cf. Table |2J) 
results in the isoscalar polarizabilities 

a s E = (13.2 ± 1.3 (stat) ± 2.1 (syst)) • 10" 4 fm 3 , 

P' M = (1.7 ± 1.6 (stat) ± 0.3 (syst)) ■ 10" 4 fm 3 , (18) 

where we assumed the same statistical errors as in Table 01 The systematic error due to the 
differing results when we use different wave functions was estimated to be around 15 %. As 
these results are in good agreement with the isoscalar Baldin Sum Rule, cf. Eq. (fT3*|) . we 
also used the central sum-rule value as additional fit constraint, obtaining 

a% = (13.1 ± 0.7 (stat) ± 2.0 (syst)) • 10~ 4 fm 3 , 

j3' M = (1.5 T 0.7 (stat) ± 0.2 (syst)) • 10" 4 fm 3 . (19) 

Motivated by the statistical imbalance between experimental data around 94.2 MeV 
and 69 MeV, we reduced in a second fit the statistics at 69 MeV, replacing the nine data 
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Figure 11: C(e 3 ) SSE (solid) and O(f) HB^PT (dashed) results with a%, f3 s M from TableH 
using the chiral NNLO wave function [2E|. The upper panels correspond to a fit of both 
polarizabilities, in the lower panels the Baldin Sum Rule (cf. Eq. (|13|)) is used as additional 
fit constraint. The grey bands are derived from our (statistical) errors. The two data points 
plotted as boxes are the effective data that we use as representative pseudo-data points for 
the data from ^UJ. The full set of data points can be seen in Fig. ® 



points given in ^U] by two representative pseudo-data points, leading to an equal weighting 
between the two energy sets. The bias-corrected fitting procedure confirms our findings 
of small values for f} s M , implying small isovector components. Our results for the isoscalar 
polarizabilities that we derive from the fit including only the two representatives of the data 
from PU] are 

a s E = (12.6 ± 1.4 (stat) ± 1.9 (syst)) • 10~ 4 fm 3 , 

P S M = (2.3 ± 1.7 (stat) ± 0.3 (syst)) • 10" 4 fm 3 . (20) 
Including the Baldin constraint we get 

a% = (12.4 ± 0.8 (stat) ± 1.9 (syst)) • 10" 4 fm 3 , 

P* M = (2.1 =f 0.8 (stat) ± 0.3 (syst)) • 10" 4 fm 3 . (21) 
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We consider these "unbiased" results, Eqs. (|2U|) and (|21|). to be the more reliable values, 
since in a straightforward fit to the existing 'yd data there is an obvious imbalance between 
the number of points at 69 and 94.2 MeV. The data base should be enlarged at higher 
energies so that unbiased fit results can be obtained. If further experiments, as planned 
at TUNL/HI7S or at MAXlab, provide additional data at energies between 70 MeV and 
the pion mass over the whole angular range, an unbiased fitting routine including data over 
this entire energy region will be possible. However, we caution that we found a very strong 
energy dependence of the 'yd cross sections in the forward direction. Therefore, we would 
recommend that any future data taken over a range of photon energies be analyzed using a 
model which incorporates this rapid energy dependence. 

A previous analysis of 'yp data within the same chiral Effective Field Theory used here 
yielded the values a p E = (11.04 ± 1.36) • 10~ 4 fm 3 , (3 P M = (2.76 =f 1.36) • 10" 4 fm 3 for the 
proton polarizabilities [2], cf. Eq. (J2J), which are consistent with experimental values [3]. 
Combining the numbers of Eq. (|2*Uj) with these results, we obtain a consistent Effective Field 
Theory determination of the neutron polarizabilities with a precision comparable to [I]: 



Eq. (|22|) does not include the Baldin Sum Rule, whereas the one-parameter fit using the 
Baldin constraint gives 



It is clear from Eqs. (j2J), (}2*2*j) . ()23|) that the isovectorial components - i.e. the differences 
between proton and neutron polarizabilities - are rather small. This finding is in good 
agreement with [I], where quasi-elastic Compton scattering off the proton and neutron was 
measured. Eqs. (j22J) and ()23j) prove that small isovectorial nucleon polarizabilities are not 
in contradiction with elastic deuteron Compton scattering data. We conclude that both 
the quasi-elastic and the elastic deuteron Compton experiments are consistent with small 
isovectorial polarizabilities. 

Furthermore we used the 0(q 3 ) HB^PT amplitudes for analogous fits, finding similar 
values for a# but larger ones for 0m, which is not surprising, as the dynamics of the resonant 
Compton multipoles is not well captured in third-order HBxPT. Therefore, the static value 
becomes large, since it must correct for the missing A resonance, leading 0(q 3 ) HB^PT to 
a disagreement with the single-nucleon Compton multipoles extracted in theories with an 
explicit A(1232), e.g. Ref. 0. Obviously, -yd scattering alone is not sufficient to investigate 
the relevant low-energy degrees of freedom in nuclear Compton scattering, but one has to 
combine information from 'yd and 'yp scattering and analyze both in the same framework. 

Finally, in future work one needs to address the issue of the sensitivity to the wave 
function (cf. Sect. I3.3|) . as one would expect that for photon energies below 100 MeV any 
effects of short distance, i.e. high-energy physics, should be able to be encoded in counter- 
terms within a well-understood power counting scheme. 




(14.2 ± 2.0 (stat) ± 1.9 (syst)) • 10" 4 fm 3 
(1.8 ± 2.2 (stat) ± 0.3 (syst)) • 10~ 4 fm 3 



(22) 



a n E = (13.8 ± 1.6 (stat) ± 2.1 (syst)) ■ 10~ 4 fm 3 , 
/3» = (1.4 t 1-6 (stat) ± 0.2 (syst)) • 10" 4 fm 3 . 



(23) 
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